Initial correlations in nonequilibrium Falicov-Kimball model 



Minh-Tien Iran 

Asia Pacific Center for Theoretical Physics, POSTECH, Pohang, Republic of Korea and 
Institute of Physics and Electronics, Vietnamese Academy of Science and Technology, 10 Dao Tan, Hanoi, Vietnam. 

The Keldysh boundary problem in a nonequilibrium Falicov-Kimball model in infinite dimensions 
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is turned on. The Kadanoff-Baym- Wagner equations for the nonequilibrium Green functions are 
derived, and numerically solved. The contributions of initial correlations are studied by monitoring 
the system evolution. It is found that the initial correlations are essential for establishing full 
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equilibrium. By examining the contributions of the initial correlations to the electric current and 
the double occupation, we find that the contributions are small in relation to the total value of those 
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strong. The neglect of initial correlations may cause artifacts in the nonequilibrium properties of 
the system, especially in the strong interaction case. 
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I. INTRODUCTION 

The theoretical description of the physical properties of 
nonequilibrium correlated electron systems is an impor- 
tant problem in condensed matter physics. In transport 
processes often correlated electron systems are driven out 
of equilibrium by switching on external fields. The sys- 
tems can be also out of equilibrium by suddenly chang- 
ing their parameters. Nonequilibrium correlated electron 
systems, which can be realized in many experiments, may 
have unusual and interesting properties. One such sys- 
tem is the quantum dot attached to two leads through 
two tunnel junctionsii^ The conductance of the dot re- 
veals a nonequilibrium Kondo effect. Other examples 
are the effects of electron correlations on the nonlinear 
current-voltage characteristics^^ Recently, experiments 
with ultracold atomic gases have made it possible to pre- 
pare initial state to a rapid change of system parameters, 
and observed remarkable subsequent dynamics as a col- 
lapse and revival of the initial phase. 6,7 

The many-body formalism for nonequilibrium sys- 
tems was developed by many people including Kubo, 8 
Schwingerj^ Kadanoff and Baying Keldyshii (see also 
Ref. [r2lfl3l for references). In particular, Kadanoff 
and Baym constructed a system of equations for the 
nonequilibrium Green functions*^ Parallel to this de- 
velopment, Keldysh also derived a perturbation theory 
for the nonequilibrium Green functions Like the Feyn- 
man perturbation theory for equilibrium systems^ the 
Keldysh nonequilibrium perturbation theory is based 
on the assumption of an adiabatic switching on of the 
many-body interactions. The assumption is necessary 
for the application of the Wick theorem, which re- 
quires a quadratic form of the system Hamiltonian at 
the initial preparation of the system. While the as- 
sumption is exactly proved in scattering theory,— its 
application to nonequilibrium many-body systems im- 
poses restrictions ! 15 ' 16 ' 17 It turns out that the assump- 



tion corresponds to the neglect of the so-called initial 
correlations i 15 ' 16 ! 17 Despite the neglect of the initial cor- 
relations the Keldysh theory is widely used in the study 
of nonequilibrium systems. Wagner unified the Feyn- 
man, Matsubara and Keldysh perturbation theories into 
a single many-body formalism in which neither a special 
form of the Hamiltonian at the initial time nor subse- 
quent time development of the system are restricted^ 7 - 
He introduced a matrix representation for the contour- 
ordered Green function and derived the Kadanoff-Baym 
equations for the nonequilibrium Green functions. In 
the Kadanoff-Baym- Wagner formalism the initial correla- 
tions are fully taken into account. While the nonequilib- 
rium formalism is well established, the role of the initial 
correlations is less attended, especially for nonequilib- 
rium strongly correlated electron systems. In particular, 
it would be desirable to test whether the initial correla- 
tions which are neglected in the Keldysh formalism are 
negligible or not. The difficulties arising in the study 
of the initial correlations are mostly due to the lack of 
getting the exact solutions of nonequilibrium correlated 
electron systems. 

In the last decade, the dynamical mean-field theory 
(DMFT) was developed] 18 ' 19 In the equilibrium case 
the theory is widely and successfully applied to study 
strongly correlated electron systems. The DMFT gives 
the exact solutions in infinite dimensions. Recently, a 
version of the DMFT for nonequilibrium systems was 
developed 31 The nonequilibrium dynamical mean- field 
theory (NEDMFT) is formally formulated on the same 
basis as of the equilibrium DMFT. Like the equilibrium 
case, in infinite dimensions the self energy of nonequi- 
librium systems becomes a local function in space. As 
a consequence, it can be self-consistently determined by 
mapping the lattice problem onto an effective problem 
of a single site embedded in a self-consistent effective 
medium. When the self-consistent equations are solved, 
the nonequilibrium Green functions are obtained and var- 
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ious physical quantities can be calculated. 

The aim of the present paper is twofold. First, we 
study the contributions of the initial correlations in a 
nonequilibrium correlated electron system. The ini- 
tial correlations are studied within truncated and self- 
consistent perturbation theories as well as within the 
NEDMFT. The Keldysh perturbation theory usually ar- 
gues for the neglect of the initial correlations. However, 
in the present paper the results obtained within the trun- 
cated and self-consistent perturbation theories show that 
the initial correlations are always finite even when the 
initial time is in the remote past limit. In the infinite 
dimension limit the initial correlations can also be ob- 
tained exactly since in this limit the NEDMFT gives the 
exact solutions. In the such way one can find under what 
circumstance the Keldysh formalism is safely applied to 
nonequilibrium correlated electron systems. The second 
aim of the present paper is to derive the Kadanoff-Baym- 
Wagner equations for the nonequilibrium Green functions 
within the NEDMFT. These equations are an alterna- 
tive to the original NEDMFT equations for the contour- 
ordered Green function^ Within the Kadanoff-Baym- 
Wagner formalism the nonequilibrium Green functions 
clearly satisfy their boundary conditions. The Kadanoff- 
Baym- Wagner formalism already includes the Keldysh 
formalism as its part, and it is suitable to study the 
initial correlations. In this paper we will examine the 
initial correlations within a Keldysh boundary problem. 
The problem works in a system which is first started 
in equilibrium and then is driven out of equilibrium by 
turning on of an external field. The model which we 
adopt to describe the system is a nonequilibrium Falicov- 
Kimball model. The Falicov-Kimball model (FKM) was 
first introduced for modelling a metal-insulator transition 
in equilibrium^ The model is one of the simplest mod- 
els for strongly correlated electron systems. The FKM 
describes conduction electrons interacting via a repul- 
sive contact potential with localized electrons. It can be 
viewed as a simplified Hubbard model where electrons 
with down spin are frozen and do not hop. Much progress 
has been made on solving this model in both exact and 
approximation ways, where all properties of the conduc- 
tion electrons in equilibrium are well knowni 22 i 23 i 24 In 
equilibrium the FKM describes a metal-insulator tran- 
sition for the homogeneous phase. 22 i 23 ' 24 The Coulomb 
interaction is divided into two ranges: the weak inter- 
action range, when the interaction strength is smaller 
than the half bare bandwidth, and strong interaction one 
otherwise. For weak interactions the system is metallic, 
and for strong ones the system is insulator. The system 
is driven out of equilibrium by a constant electric field. 
The electric field is switched on at a some time after 
the initial preparation of the system in equilibrium. This 
nonequilibrium FKM was introduced by Freericks et al in 
the study of the Bloch oscillations in the electric current 
within the NEDMFT^ In the present paper we derive 
the Kadanoff-Baym- Wagner equations of the NEDMFT 
for the nonequilibrium FKM, and calculate the contri- 



butions of the initial correlations to the electric current 
and the double occupation. It is found that the contri- 
butions of the initial correlations to the electric current 
and the double occupation are small in relation to the full 
value of those physical quantities in the weak interaction 
case, and significantly increases in the strong interaction 
case. However, without the initial correlations the sys- 
tem cannot restore full electron correlations even before 
the turning on of the electric field. The neglect of the 
initial correlations may cause artifacts in the nonequilib- 
rium properties of the system. 

The paper is organized as follows. In Sec. II we present 
the Kadanoff-Baym- Wagner nonequilibrium formalism, 
and describe the nonequilibrium FKM. In the next three 
sections we present the studies of the model within the 
truncated and self-consistent perturbation theories, and 
the NEDMFT. The last section is the conclusion. 



II. FORMALISM AND MODEL 

We consider the Keldysh boundary problem in a 
nonequilibrium system which is first prepared in equilib- 
rium, and then is driven out of equilibrium by switching 
on of an external field or by a sudden change of its pa- 
rameters. Specifically, at an initial time to the system 
is prepared in equilibrium which is defined by the equi- 
librium Hamiltonian H cq and temperature 1//3, and at 
time t = (to < 0) an external field is switched on or its 
parameters are suddenly changed. Usually, in nonequi- 
librium systems the time translational invariance is not 
valid, and the Green functions, which are employed for 
studying the physical properties of the systems, depend 
on the two time variables. The nonequilibrium formalism 
works with the so-called contour-ordered Green function 
which is defined for the time variables on the Kadanoff- 
Baym contouri 10 ' 11 ! 12 ' 13 The Kadanoff-Baym contour is 
shown on Fig. [T] The contour starts in the initial time 
to, runs out to maximal time t m , then returns to the 
initial time, and finally moves parallel to the negative 
imaginary axis a distance /3. At the initial time to the 
system is always in equilibrium. The Kadanoff-Baym 
contour is suitable for deriving the Dyson equation for 
the contour-ordered Green function. Keldysh also intro- 
duced a similar contour for the contour-ordered Green 
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FIG. 1: Kadanoff-Baym contour for the two-time Green func- 
tions in nonequilibrium. 
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function^ The Keldysh contour is basically the same as 
the Kadanoff-Baym contour, but it neglects the last con- 
tour branch parallel to the imaginary axis and limits to 
to minus infinity. The neglect of the last branch of the 
contour corresponds to the neglect of initial correlations. 
The contour-ordered Green function is defined by 



G c {i,j\m = -i(T cCi (t)c](F)) 



+ib 



(1) 



where cj (cj) are the creation (annihilation) operators for 
electrons at site i. The time evolution of the operators 
on the Kadanoff-Baym contour is defined in the Heisen- 
berg picture. T c is the time ordering on the Kadanoff- 
Baym contour and it is defined via the contour step func- 
tion 8 c (t,t')- 6 c (t>t') equals to 1 if i lies after P on the 
contour, and it equals to otherwise. The averages in 
Eq. (JTJ) are the statistical average over the equilibrium 
Hamiltonian H oq at temperature l/f3. The Kadanoff- 
Baym contour consists of three time branches: the first 
branch is chronological, the second one is antichronolog- 
ical, and the last one is parallel to the imaginary axis. 
Thus, we can represent the contour-ordered Green func- 
tion by a 3 x 3 matrix G a7 (i, F), where i locates on the 
a-th branch, and F locates on the 7-th branch. In the 
such way, the contour-ordered Green function has nine 
component Green functions, however they are not in- 
dependent. Wagner reduced the matrix representation 
of the contour-ordered Green function to a matrix form 
of six component Green functions, and five of them are 
independent^ The Wagner matrix representation for 
the contour-ordered Green function can be written as 
follows^ 



G R G K 
(V j G A 
V2G^ G M 



(2) 



where 



G R (t,F) 
G A (t,t') 
G K (t,t') 
G M (t, t') 
G\tS) 



G n (t,t')-G 12 (t,t') 

-ie{t-F){{ Ci {t),c){F)}), 

G 12 (t,t')-G 22 (t,t') 

ie{F-t)({ Ci {t),c]{F)}), 

G 12 (t,t') + G 21 (t,t') 
-i([ci(jt),c]W]), 

G 33 (t -iT,t -lT') 

-i(T T Ci(t - ir)c](t - ir')), 
G 13 (t,t -ir') 
^<c 3 t (^ -^r , )c ^ (^)), 
G 31 (t Q -iT,t') 

-i(ci{to-iT)c](F)), 



where t, F are real times, and < r, t' < (3. In the 
above equations we have used the commutator symbol 
[A, B] = AB — BA, and the anticommutator symbol 
{A, B} = AB + BA. G R{ - A) are the retarded (advanced) 
Green function, and G K is the Keldysh Green function. 
These Green functions are defined totally on the real 
time axis. G M is the Matsubara Green function and is 
defined on the imaginary time branch of the Kadanoff- 
Baym contour. Note that the above definition of the 
Matsubara Green function differs from the standard one 
by factor The Green functions G^ and G^ have one 
time variable on the real time axis, and the other vari- 
able on the imaginary time branch. 2 - They do not have 
a specific name, however we will refer them to the right 
and left time mixing Green function, respectively. The 
left corner 2x2 matrix in the Wagner matrix represen- 
tation in Eq. ([2]) is the Keldysh representation of the 
nonequilibrium Green functions in the Keldysh perturba- 
tion theory^ The Matsubara Green function is just the 
equilibrium Green function at temperature 1/(3. It cou- 
ples with the Keldysh Green function through the time 
mixing Green functions. If the time mixing Green func- 
tions are neglected the Kadanoff-Baym- Wagner formal- 
ism reduces to the Keldysh formalism. 

Within the Wagner matrix representation the Dyson 
equation for the nonequilibrium Green functions can be 
written in the standard form like in the equilibrium case 



G = Go + G • S • G, 



(3) 



where Go is the bare Green function, and £ is the self 
energy. The self energy is also written in the Wagner 
matrix representation 



^=|0 £ A 

\/2£r E M 



(4) 



Note that in the Dyson equation the product symbol 
• denotes not only the matrix multiplication, but also 
the integration over the time variables. We also omitted 
other variable notations such as of momentum or spin to 
simplify the equation writing. We introduce the inverse 

matrix Green function G by the standard definition 



G.G= I. 



(5) 



The inverse matrix Green function is also presented in the 
Wagner matrix representation. One can find its elements 
by explicitly writing the component equations of Eq. (|5|) 



G ■ G 



K 



G R ' A ■ G R ' A = i, 
G M * G M = I, 
G K ■G A + 2G ] *G r = 0, 



G R ■ G 1 + G 1 * G 
& ■ G 1 



0. 



G M *G^ 



(6) 
(7) 
(8) 
(9) 
(10) 
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Here the dot and star products are the integrations over 
the real time and the imaginary time variables, respec- 
tively, i.e., 

(A ■ B)(i,F) = [ m dhAfct^Bih,?), 

J to 

r 

(A*B)(t,t') = -i driA(t,ri)B(ri,t'). 
Jo 

The symbol 1 is just the delta function of the time vari- 
ables. For real time variables it is 8{t — t'), and for imagi- 
nary time variables it is iS(r—T'). One can view G R I A I M 
as the inverse matrices of G R ' A I M in continuous time 
variables. However, G K /~\/t are not inverses of the corre- 
sponding Green functions. The Dyson equation ([3]) can 
be rewritten as follows 



Go • G = 1 + E • G, 



(11) 



where Go is the inverse matrix of Go , and its elements can 
be found from Eqs. (f6|)- (|10p for the bare Green functions. 
Equation (llip can be written in the explicit form for the 
component Green functions 



QR/A , g r /A 

Gq 1 -k G M 
G R • Gl 



= i + j: r / a -g r / a , 



= 1+E M *G 



(12) 
(13) 



G n M 



G R -G 



K 



= Y, R -G^ -G )*G m , (14) 

= £ M * G^ + (E^ — Gq) • G A , (15) 



E • G 



K 



Yi K — Gn ) ■ G 



K \ 



+2(£l -G\) *G r . 



(16) 



Equations (ri"2"]) - (fT6")) are just the Kadanoff-Baym equa- 
tions for the nonequilibrium Green functions written in 
the Wagner representation. In the standard Kadanoff- 
Baym equation o 10 ! 17 the inverse bare Green functions 
are written in the form of differential operators, and 
these differential equations also require additional bound- 
ary conditions. In the Kadanoff-Baym- Wagner equations 
(fT2"]) - (|16p the inverse bare Green functions have their ex- 
plicit forms and they already satisfy their boundary con- 
ditions. Instead of differential-integral equations in the 
Kadanoff-Baym formalism of the contour-ordered Green 
function, the Kadanoff-Baym- Wagner equations (fT5|) - 
([11]) are just only the integral equations. Once the self 
energy is computable the Kadanoff-Baym- Wagner equa- 
tions (fT2"|) - (fT6")) can be solved. First we solve Eqs. (TT2"|) - 
(TlU]) for the retarded, advanced and Matsubara Green 
functions. These equations can be solved independently. 
Certainly, the advanced Green function can be quickly 
obtained from the retarded Green function by the rela- 
tion 



G A (t,t')=[G R (t',t)]' 



(17) 



Moreover, the Matsubara equation (flU]) is the equilibrium 
equation and we can also use the equilibrium techniques 



to calculate the Matsubara Green function. Next we use 
the retarded, advanced and Matsubara Green functions 
as the inputs and solve the next two equations for the 
time mixing Green functions. Finally, we solve the last 
equation for the Keldysh Green function. The Kadanoff- 
Baym- Wagner equation for the Keldysh Green function 
([in]) can be rewritten as 



G 



K 



(1 + G R ■ T, R ) ■ Gq ■ (i + T, A ■ G 
+G R ■ Yi K ■ G A — 2G R 



Gn * Gn^ -k Gn 



-(El -G )*G M *(l]r-G ) ■ G A . (18) 

Here we have used Eqs. ©, (JTOJ for Gjf , Gq and 
Eqs. (TT2"1) . (Ti"5)) for the retarded (advanced) and time mix- 
ing Green functions. If the time mixing Green functions 
are neglected, the Kadanoff-Baym- Wagner equation (JTSJ) 
is reduced to the Keldysh equation 



G K = (1 



G R ■ E fl ) • Gq ■ (i 

K 



G A ) 



-G ■ E • G 



(19) 



The Keldysh formalism neglects the contributions gener- 
ated from the dynamics of the system in the imaginary 
time branch of the Kadanoff-Baym contour. Since at the 
initial time to the system is prepared in equilibrium with 
full interaction, the neglected contributions are correla- 
tions of electrons between the initial time and an ad- 
vanced time. Indeed, if we neglect the correlation effects 
of the Matsubara and the time mixing Green functions 
(i.e., E M /1/T = 0), the last term in Eq. ((T5J) vanishes, and 
we again obtain the Keldysh equation. The neglected 
contributions are called initial correlations! 13 ' 15 i 16 ' 17 The 
initial correlations distinguish between the Kadanoff- 
Baym- Wagner and the Keldysh formalisms. One can no- 
tice that the equations for the retarded and advanced 
Green functions are decoupled from the system of equa- 
tions, hence the nonequilibrium density of states remains 
the same in both the Kadanoff-Baym- Wagner and the 
Keldysh formalisms. The initial correlations do not affect 
the nonequilibrium density of states. They affect only 
the nonequilibrium distribution function. Thus the ini- 
tial correlations give contributions only to physical quan- 
tities which depend on the nonequilibrium distribution 
function. 

The model we will study is the FKM with external 
electric field turned on at t — 0. At the initial time to 
the system is prepared in equilibrium with temperature 
1/(3 and the FKM Hamiltonian 

i,j i i 

+Uj24*fifi, (20) 

i 

where c] (cj) are the creation (annihilation) operators 
for conduction electrons at site i, and /■ (/j) are the cre- 
ation (annihilation) operators for localized electrons at 
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site i. Jij is the hopping matrix of conduction electrons, 
and it is equal to J for nearest neighbor sites, and it is 
otherwise. U is the strength of the interaction between 
the conduction and localized electrons. \x and Ef are the 
chemical potentials of the conduction and localized elec- 
trons, respectively. In this paper we will only consider 
the half filling case. It turns out that in the half filing 
case fi = —Ef = U/2. At time t = a spatially uniform 
electric field is turned on. We choose the gauge with van- 
ishing of the scalar potential for the electric field. As a 
result the electric field is described by a spatially uni- 
form vector potential A(t) = —9(t)Et. The electric field 
couples to the conduction electrons through the Peierls 
substitution for the hopping matrix 



Ji 



Jij exp 



i a I drA(r, t) 



Jij exp [ — zeA(t)(Rj — 



(21) 



By replacing the hopping matrix in Hamiltonian in 
Eq. (f2T)|) by Eq. (f2"Tj) we obtain full nonequilibrium 
Hamiltonian of the system. This nonequilibrium FKM 
was introduced by Freericks et al in the study of the 
NEDMFT^S, The considered nonequilibrium FKM dif- 
fers from the equilibrium FKM only by the bare energy 
spectra 

d 

e{k,t) = e(k-eA(i)) = -2J^cos(fc, - eAj(*)),(22) 

i=l 

where d is the space dimension of the system. We will 
consider the case when the electric field lies along the 
elementary cell diagonal 

A(t)=A(t)(l,l,...,l). 

In this case the bare energy spectra becomes 

e(k, t) = cos(eA(t))e(k) + sin(eA(i))e(k), (23) 

where 

d 

e(k) = -2J^cos(fc,), 

i=i 
d 

e(k) = -2J^sin(fci). 



In the limit of infinite dimensions d — > oo the bare density 
of states has a double Gaussian form 



p{e,e) = po(e)po(e), 



(24) 



where po(s) = exp(— £ 2 )/y / 7t. Here we have used J* = 
J\[d as the unit of energy. 

The nonequilibrium bare Green functions can be found 
from the equations of motion. The equation of motion 
for the retarded Green function reads 

[idt + mo - e(k, t)] G*(k\t, t') = 6(t - t'), 



where (1q is the chemical potential of the noninteracting 
conduction electrons. At half filling (j,q = 0. With the 
boundary condition G?j^(k|t, t) = —i we can find the bare 
nonequilibrium retarded Green function 

G$(k\t,t') = -i6(t-t')e i ^ t - t 'h- i ft' dtl<Ktl) . (25) 

Similarly, one can find the bare advanced, Keldysh and 
Matsubara Green functions 

G$(k\t,t') = i6(t' -t)e i ^ t - t ')e- i ti dtis(k ' tl \(26) 
G«(k\t,t') = *[2/( £ (k)-Mo)-l]e^ (t -*' ) 

e -i/ t **i e (k,ti)^ ^ 

G M (k|r,r') = -i[e(r-r')-f(s(k)-n )] 



3 -(e(k)- Aln )(r-r') 



(28) 



where /(e) = l/(exp(/3e) + l) is the Fermi-Dirac distribu- 
tion function. The bare right time mixing Green function 
can be found from the equation of motion 

[id t +fx -e(k,t)]Gl(k\t,T') =0, 

with the boundary condition Go(k|t ,T') = Go J (0,t'). 
We obtain 



G](k\t, T ') = -i[0(- T ')-/( £ (k)- Mo )]e 
= iG^(k%to)G^(k\0,T'), 



(e(k)— Ato)r' 



(29) 



since t > to. Similarly, the bare left time mixing Green 
function is 

Gi(k\r,t') = -i[6( T ) - f(e(k) - v )]e-WV-^ 

e vo(t'-t ) e i r mmkm) 

= -tG" o I (k\T,0)G^(k\t ,t'). (30) 



The nonequilibrium bare Green functions clearly satisfy 
their boundary conditions. When the self energy is com- 
putable, it together with the bare Green functions fully 
determine the nonequilibrium Green functions via the 
Kadanoff-Baym- Wagner equations (P~2|) - fTBj) . In the next 
sections we will solve the Kadanoff-Baym- Wagner equa- 
tions with the self energy calculated within the truncated 
and self-consistent perturbation theories as well as within 
the NEDMFT. 



III. TRUNCATED PERTURBATION THEORY 

In this section we calculate the nonequilibrium Green 
functions and the electric current within the truncated 
perturbation theory of second order in U. The pertur- 
bation calculations were previously performed within the 
Keldysh nonequilibrium perturbation theory^ where the 
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initial correlations are neglected. The purpose of this sec- 
tion is to find the contributions of the initial correlations 
to the electric current within the truncated perturbation 
theory. 

In the half filling case the first-order perturbation con- 
tributions to the self energy vanish^ The second-order 
self energy can be found by expanding the contour- 
ordered Green function to second order in U . One can 
find 2 ^ 

££(i, = U 2 n f (l - n f )j- ]T G?(k|t, t>), (31) 

k 

where a = R, A, M, K,~\,\, and n/ = 1/2 is the density of 
the localized electrons at half filling. Within the second- 
order perturbation the self energy does not depend on 
momentum. This feature is similar to the DMFT where 
the self energy is a function of time variables only. Using 
the Dyson equation ([3]) we can obtain the nonequilibrium 
Green functions up to second-order in U 

G%(k\t,t')=G *(k\t,t , ) + AG%(k\t,t'), (32) 

where 

AG* /A = G* /A ■ ■ G« /A , (33) 

AGf = G^*£f *G^, (34) 

A/^l" — fiA I fiM j. v< T riA 

IAU 2 — U Q ■ 2j 2 • U -t- U Q *2j 2 -U 

+G M *Ef*4, (35) 

Alj 2 — °0 ^2 ' "0 + U ' Zj 2* u 

+Gj*£ 2 M *G^, (36) 
a s~tK _ r<R r<K x_nR r'- 4 

+Gq ■ Yi A ■ G A 4- 2G^ • £ 2 * Gq 
+2Gj*S|-G^ + 2Gj*Ef *G£. (37) 

Here, for simplicity we omitted the variable notations in 
the Green functions and the self energy. In comparison 
to the Keldysh perturbation theory, the Kadanoff-Baym- 
Wagner perturbation expansions of the retarded and ad- 
vanced Green functions remain unchanged^ However, 
the second-order perturbation expansion of the Keldysh 
Green function is different. It consists of two parts. The 
first part is the first three terms in Eq. ([37]) which are also 



the perturbation contributions within the Keldysh per- 
turbation theory, 2 ® and the second part is the remaining 
last three terms which are additional contributions gen- 
erated from the initial correlations. The Keldysh pertur- 
bation theory neglects the second part. 

The electric current can be calculated by evaluating 

k 

where «i(k) = J*sm(ki)/\/d is the velocity component, 
and G < (k|t,t) is the equal time lesser Green function, 
which can be calculated from the Keldysh Green function 
by the relation 

G<(k|M) = i(G A '(k|M) + z). (39) 

When the electric field lies along the diagonal, all compo- 
nents of the electric current are equal, and the magnitude 
of the current is 

I{t) = Vdli{t). (40) 

By inserting the second order perturbation expansions of 
the self energy in Eq. (|3"Tj) and of the Green functions in 
Eq. ((221) into the current formulas in Eqs. ([55 )1 - ((ID )) , we 
obtain the electric current up to order U 2 

h(t) = I (f) + AI 2 (t), (41) 

where Io(t) and A/ 2 (i) is the zeroth and second order 
contributions to the current. The zeroth order current is 

I (t)=j J dep (e)ef(e)sm(eA(t)), (42) 

where jo = e/yd. It is the electric current in the non- 
interaction case. It exhibits the Bloch oscillations with 
period 2ir/E, and its amplitude is independent on time. 
In the noninteraction case the Bloch oscillations of the 
current occur when the noninteracting electrons move in 
a lattice under a constant electric field. In this case the 
system is a perfect conductor, the periodicity of the lat- 
tice restricts the wave vector to lie in the first Brillouin 
zone that leads to the oscillations of the current. 



After some analytical calculations we also obtain the second-order perturbation contributions to the current strictly 
in the half filling case 



AI 2 {t) 
A/f(t) 



A/f(i)+A/| c (t), 

jj2 rt rti 
T 



(43) 



to 



jo~r I dti I dt 2 / dep(e)tauh( — ) exp 



:cos (eC(t2,h)) sin(eA(t)) 



■ sin e 



C(t2,tl))S(t 2 ,h) cos(eA(t)) 
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-Jo 



U 



2 rt 



16 



dtx j dt 2 J dep(e) tanh ex P 



sin ( 



(e<7(*a, *i)) C(t 2 ,h) sm(eA(t)) - S(t a , *i) cos(eA(t)) 



3o- 



U 2 



s' - e 
U 2 



S(t 1: t ) cos(eA(t)) + i(e - e') sm(eA(t)) 



+j — JdeJ ds'p(e)p(e')f(e)f(~e)esm(eA(t)) 
/(£ ) — 



(e-e>) 



A2 



(5-5') 



M2 



(44) 



dti / rfe / rf£'p(e)|o(e') exp [ - -S 2 (ii, t ) + «(e - e')C(*i, *o) 



(45) 



r 



Here in order to simplify the expression, we have in- 
troduced the functions 

C(t 2 ,h) = / dt' coa(eA(t')), 
Jti 

S(t 2 ,h) = [ 2 dt'sm(eA(t')). 
Jti 

Like the Keldysh Green function the second-order current 
also consists of two parts, AI 2 (t) and AI 2 c (t). The first 
part AI 2 (t) is the second-order contributions within the 
Keldysh perturbation theory^. The second part AI 2 c (t) 
is the contributions of the initial correlations. The second 
part is beyond the Keldysh perturbation theory. 

In Fig. [2] we plot the electric current calculated within 
the Kadanoff-Baym- Wagner and the Keldysh perturba- 
tion theories up to second order in U. For comparison 
we also plot the exact result which is obtained by per- 
forming the NEDMFT calculations (see Sec. V). It shows 
that the current has the Bloch oscillations with period of 
2ir/E like the current in the noninteraction case^S How- 
ever, the amplitude of the current varies with time. The 
perturbation theories give reasonable results for times 
smaller than ~ 2/U. The Kadanoff-Baym- Wagner per- 
turbation theory overestimates the current, while the 
Keldysh perturbation theory underestimates it. Figure^ 
also shows that the Kadanoff-Baym- Wagner perturbation 
result is closer to the exact solution than the Keldysh 
perturbation one at times right before the perturbation 
theories are broken down. In Fig. Owe also plot the ini- 
tial correlation contribution AI 2 c (t) to the current. This 
part of the current also oscillates with the same period 
as of the full current. In Fig. [3] we plot the initial cor- 
relation part of the current for various values of U and 
E. In contrast to the full current, the amplitude of the 
initial correlation part does not significantly vary with 
time. Since the initial correlation contribution is cal- 
culated within the second-order perturbation theory, its 
amplitude is proportional to U 2 , and almost independent 
on the electric field. Usually, the Keldysh perturbation 
theory argues that the initial correlations vanish when 
the initial time approaches to minus infinity. In Fig. 0] 



we plot the current and its initial correlation part at a 
fixed time as functions of the initial time to. It shows that 
both the current and its initial correlation part quickly 
approach to constant values when |£q| increases. Even 
for to — — 5 the current and its initial correlation part 
already reach the constant values. The initial correla- 
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FIG. 2: (Color online) The time dependence of the electric 
current calculated within the Kadanoff-Baym- Wagner (KBW, 
red line) and the Keldysh (K, blue line) perturbation theory 
(PT). The exact NEDMFT calculation result and the initial 
correlation contribution to the current are presented by the 
black and green lines, respectively. The model parameters 
U = 0.5, (3 = 10, t = -10, and E = 1 (E = 2) for upper 
(lower) panel. 



8 




10 



15 



FIG. 3: (Color online)The second-order initial correlation 
contribution to the current as a function of time for various 
U and E (to = -10, (3 = 10). 



tions never vanish, even when t a — > — oo. Thus the 
Keldysh perturbation theory always neglects the nonvan- 
ishing initial correlations. However, within the truncated 
perturbation theory both the Keldysh and the Kadanoff- 
Baym- Wagner formalisms only qualitatively describe the 
physical properties when the perturbation theory works. 
The initial correlations do not qualitatively change the 
perturbation results. Thus, the use of the Keldysh per- 
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FIG. 4: The electric current (solid line) calculated within 
the Kadanoff-Baym- Wagner perturbation theory (KBW PT) 
and its initial correlation part (dotted line) at time t = 2 as 
functions of the initial time to for U = 0.5, E = 1, and /3 = 10. 



turbation theory is still convenient in the nonequilibrium 
study due to its simple system of equations. 



IV. SELF-CONSISTENT PERTURBATION 
THEORY 

In this section we perform the self-consistent perturba- 
tion calculations for the electric current. Instead of the 
standard perturbation calculation in Eq. (j3ip . we take a 
self-consistent approach by dressing the bare Green func- 
tions in the calculation of the self-energy, i.e. 

E a (t, f ) = U 2 n f (l ~n f )±Yl GQ ( k l*> 0- (46) 

k 

In this approximation only the Green functions of the 
conduction electrons are dressed. The Green functions of 
the localized electrons are kept local, thus their contri- 
butions to the self energy of the conduction electrons are 
just — 7t/). We solve the Kadanoff-Baym- Wagner 
equations p2p - (fH)l) with the self energy determined by 
Eq. pp| . In order to solve these equations we adopt 
the discretization method which was employed by Fre- 
ericks et al in solving the NEDMFT equations32! We 
discretize the time variables with step At for real time 
t and At for imaginary time r. As a result the real 
time domain is divided into L points, and the imagi- 
nary time domain (3 is divided into M points. Thus, any 
function of two time variables A(t, t') becomes a matrix 
Aij = A(ti, tj), where U = t and tj — t' . Integration over 
time can be approximated by the rectangular integration 
rule 

J i 

where At = At for real time integration, and At = —iAr 
for imaginary time integration. Thus the time integra- 
tion becomes a matrix multiplication. The inverse of the 
continuous matrix function 



dt 1 A(t,t 1 )A- 1 (t 1 ,t') = 5{t-t') 



in the discretization approach becomes 



At ^2 A{t% , ti)A~ x (ti , tj 



u *3 

At 



(47) 



Thus, in the discretization approach the Kadanoff-Baym- 
Wagner equations become the matrix equations which 
can be solved numerically. The time discretization is 
a numerical approach which approximately solve the 
Kadanoff-Baym- Wagner equations. It becomes exact 
only for At — * 0. Nevertheless, it was shown that the 
discretization approach is an efficient way to solve the 
nonequilibrium Green function equations! 20 ! 2 Note that 
the Kadanoff-Baym- Wagner equations P^)l - (fTB|) differs 
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from the contour-ordered Green function equation^ Nu- 
merically, here we have to solve the equations of matri- 
ces with size L x L, L x M, and M x M, instead of 
matrices of size (2L + M) x (2L + M) in the contour- 
ordered Green function equation. It reduces the matrix 
size and computation time. However, here we have to 
solve five equations with additional matrix multiplica- 
tions. The inverse bare Green functions Gq are cal- 
culated from Eqs. (f6))- (fl"0|) with the inputs of the bare 
Green functions in Eqs. (p25|) - (P0)) . Within the discretiza- 
tion accuracy, these inverse bare Green functions are cal- 
culated exactly. They satisfy the boundary conditions. 



For instance, the Matsubara Green function has the an- 
tiperiodic property in the time variable, or the Keldysh 
Green function satisfies Gff(t ,t ) = 2(7^(0, 0+) - i. In 
the contour-ordered Green function approach, the inverse 
bare Green function contains a time differential operator 
and it is also approximately discretized. In the present 
approach the inverse bare retarded and advanced Green 
functions are numerically calculated from their bare func- 
tions by the discretization inverse relation in Eq. (|4T[) . 
The inverse Matsubara Green function in the discretiza- 
tion form can be analytically obtained 



G^(k) 



At 2 



{ 1 

_ e -e(k)Ar 





e -e(k)Ar \ 





-e(k)Ar l q 



V 







= -e(k)Ar 



1 



(48) 



Here we have taken into account /io = at half filling. This inverse bare Matsubara Green function corresponds to the 
bare Matsubara Green function with fixed diagonal elements Gq' 1 (t, t) = —i(l — /(e(k))). It is suitable for calculating 

the left time mixing Green function because of the boundary condition Gq(t, to) = Gq 1 (t, 0) for r > 0. The right 

time mixing Green function has the boundary condition G (io, t') = (0, t') for r 1 > 0, and the Matsubara Green 
function suitable for its calculations has the diagonal elements Gq 1 ^, t) — if(e(k)). The corresponding inverse bare 
Matsubara Green function has the matrix form 



/ 1 




G^(k) 







£(k)Ar 



_ e e(k)Ar 
1 







_ e e(k)Ar 
1 





_ e e(k)Ar 



\ 



3 e(k)AT 
1 



(49) 



For a definiteness we also use this Matsubara Green func- 
tion for calculating the Keldysh Green function. The 
Green functions G$^^ can be also analytically ob- 
tained. From Eqs. (|6j)- (|10p for the bare Green functions 
in Eqs. one can show that 

Gl(t,T>) = 6(t-ta)6(r'), (50) 
G[{r,t') = -5(r)S(t'-t ), (51) 
G*(t,i') = i5(t-t )6(t' -t ). (52) 

In the Keldysh formalism, when the time mixing Green 
functions are neglected, the Green function Gq is little 
changed. One can obtain 

G^(M') = -<[2/(e(k))-l]5(t-t )J(t / -to).(53) 

In order to solve the Kadanoff-Baym- Wagner equa- 
tions ([T2 ) -([T5 ]) . first we solve Eqs. O and ([T3J) for 
the retarded (advanced) and Matsubara Green functions, 



and then we find the time mixing Green functions from 
Eqs. fT4")) - (fl"5")) . Finally, the Keldysh Green function is 
calculated from Eq. (fTtJ)) . We use iterations to solve 
each equation. When the noncquilibrium Green func- 
tions are obtained, we can compute the electric current 
by Eq. (|38p . The momentum summation in Eq. (I46p 
or ([55)1 indeed is the integration with the double Gaus- 
sian density of states in Eq. (|24]) . and we use a Gaussian 
quadrature to calculate it. Typically, we use 51 points 
for the Gaussian quadrature. In the next section we will 
discuss this type of integrations in a more detail. The cal- 
culated current converges with At well. We can obtain 
reliable results at At — > by using a Lagrange interpola- 
tion formula. Typically, we use a quadratic interpolation 
to obtain the current in the continuous limit. In Fig.[5]we 
plot the electric current obtained within the Kadanoff- 
Baym- Wagner self-consistent perturbation theory. For 
comparison we also plot the exact NEDMFT calculation 
result (see also the next section). Figure [5] shows that the 
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FIG. 5: (Color online) The time dependence of the electric 
current calculated within the Kadanoff-Baym- Wagner self- 
consistent perturbation theory (KBW SCPT, red line). The 
exact NEDMFT calculation result and the initial correlation 
contribution to the current are presented by the black and 
green lines, respectively. The current was already scaled 
with a quadratic extrapolation (At = 0.1, 0.065, 0.05, and 
At = 0.1). The model parameters U = 0.5, to = —5, /3 = 10, 
E = 1 (E — 2) for upper (lower) panel. 



self-consistent perturbation theory gives very good re- 
sults for time smaller than 2/U. In comparison with the 
truncated perturbation theory, the self-consistent pertur- 
bation theory gives reasonable results in a wide range of 
the time variable. The current obtained within the self- 
consistent perturbation theory also oscillates with time, 
and is damped to zero value. Even for large electric fields 
(for instance, E = 1) the time damping of the current is 
still observed in the self-consistent perturbation results 
like the exact solution. However, for larger electric fields 
(for instance, E — 2), the self-consistent perturbation 
theory cannot reproduce the beat behavior of the cur- 
rent, as shown in the lower panel of Fig. [5l It shows that 
the self-consistent perturbation theory may produce ar- 
tifacts, especially for noncquilibrium steady state. How- 
ever, this happens only for very strong electric fields. We 
define the initial correlation contribution to the current 
by the difference of the currents calculated within the 
Kadanoff-Baym- Wagner and the Keldysh self-consistent 
perturbation theories. In contrast to the results of the 
truncated perturbation theory, the initial correlation part 
of the current is damped with time, and its amplitude is 
significantly smaller. We plot the initial correlation con- 
tribution to the current for various values of U and E 
in Fig. [6j It shows that the amplitude of the initial cor- 
relation part is not scaled with U 2 . For long time limit 
the initial correlation contribution to the current van- 
ishes. In this case the Keldysh and the Kadanoff-Baym- 




FIG. 6: (Color online) The time dependence of the ini- 
tial correlation contribution to the current calculated within 
the Kadanoff-Baym- Wagner self-consistent perturbation the- 
ory for various U and E. The current was already scaled 
with a quadratic extrapolation (At = 0.1, 0.065, 0.05, and 
At — 0.1). The other model parameters to — —5, f3 = 10. 



Wagner formalisms give the same steady state. However, 
it may be an artifact, especially for very strong electric 
fields when the exact current exhibits the beat behavior. 
This also indicates that the self-consistent perturbation 
theory may not work well for very strong electric fields. 
Nevertheless, within the self-consistent perturbation the- 
ory the Keldysh formalism qualitatively gives almost the 
same results as of the Kadanoff-Baym- Wagner formal- 
ism. The initial correlations do not qualitatively change 
the perturbation results. 



V. NONEQUILIBRIUM DYNAMICAL 
MEAN-FIELD THEORY 



In this section we present the NEDMFT through the 
Kadanoff-Baym- Wagner representation. The NEDMFT 
was proposed by Freericks et al and it is based on the 
same idea of the DMFT in equilibrium^ The NEDMFT 
has the same principle features of the equilibrium DMFT. 
It becomes exact in the infinite dimension limit. In infi- 
nite dimensions the self energy is purely local in space. It 
can be determined by mapping the lattice problem into 
an effective problem of a single site embedded in a self- 
consistent effective medium. The effective medium can 
be represented by a Green function Q which is determined 
by the Dyson equation 



(54) 



where G = J^kGO^/N. From this equation we can 
find the components of the effective medium Green func- 
tion in the Kadanoff-Baym- Wagner representation like 
Eqs. (Hl])-([llD. We obtain 



G R/A + E H/A 



(55) 
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gM 


= G M + £ M , 


(56) 




= si -G k -GUg m , 


(57) 




= -G M *Gl -G A , 


(58) 


9 K 


= H K — G R ■ G K ■ G A 






+2(£l ~&)*G^ -G A , 


(59) 



where Q and G are the inverse matrices of Q and G, re- 
spectively. Equations (|55|) -(|59 |) are the Kadanoff-Baym- 
Wagner equations for determining the effective medium 

Green function Q. Once the effective medium Green 

function Q is determined we can compute the single-site 
Green function. In the homogeneous phase we obtain 20 

Gim P = (1 - n/)Qo + nfQi, (60) 

where nf is the localized electron density, and Qi with 
I = 0, 1 satisfy the following equation 

\§ + A/j, - IU] »Qi = 1. (61) 

Here A/i = /i — (1q. One can find explicitly the compo- 
nents of Qi by using the inverse equations (|6|l- (fT0|) 



Ft/A _ 


[Q R / A +Ag?' A ]-\ 


(62) 


Qf = 


[q m + a~ 9 ?]-\ 


(63) 


q] = 


*Qf, 


(64) 


q\ = 


-Ql**&-Qt, 


(65) 


Qf = 


~Qf -G K -Q? 






+2Qf +Q? 1 -Qf 


(66) 



where Agf = g a (no + A^i — IU) — g a (^) with a = 
R, A, M, and g a (x) is the inverse of the bare Green func- 
tion g a (x) of a pure noninteracting single site with zero 
energy level and the chemical potential x. Note that in 
the numerical calculations when we make the discretiza- 
tion of the time variable, the quantity Afi — IU does not 
lie in the diagonal of the matrices of the inverse retarded 
(advanced) or Matsubara Green functions. It lies in the 
first subdiagonal of the matrices like in Eqs. (|48p. (|49| . In 
the such way we can compute the single-site Green func- 
tion Gimp. However, it is applicable only for nonequilib- 
rium FKM. For other models such as the Hubbard model 
one may adopt different techniques to solve the effective 
single-site problem. 

The self-consistent condition requires that 

Gimp = G. (67) 

With this self-consistent condition when the effective 
single-site problem is solved we can again compute the 
self energy from the Kadanoff-Baym- Wagner equations 
(f55|) - (|59p . When the self energy is obtained the full lattice 
Green functions are calculated from the Kadanoff-Baym- 
Wagner equations (|T2]) - (P~H|) . Thus we obtain a closed sys- 
tem of equations for the nonequilibrium Green functions 



in the NEDMFT. Like the previous section, first we solve 
the set of equations of the retarded (advanced) and Mat- 
subara Green functions. Then use the obtained Green 
functions to solve the set of equations of the time mixing 
Green functions. Finally, we compute the Keldysh Green 
function from its set of equations. We use iterations for 
finding each Green function. Numerically, we employ the 
discretization method which was described in Sec. IV to 
solve the Kadanoff-Baym- Wagner NEDMFT equations. 
In equilibrium the FKM describes a metal-insulator tran- 
sition for the homogeneous phase d 1 ! 22 ' 23 The Coulomb 
interaction is divided into two ranges. For weak inter- 
actions (U < \/2) the system is metallic, and for strong 
ones (U > V2) the system is insulator. We will study 
the two cases separately. 

As a benchmark we apply the Kadanoff-Baym- Wagner 
NEDMFT to the equilibrium FKM at half filling. The 
DMFT results of the FKM at equilibrium can be also 
obtained by solving the DMFT equations in frequency* 2 ^ 
As noted in Sec. IV, the summation over momentum is 
replaced by integration with the double Gaussian density 
of states in Eq. (|24p . and we use a Gaussian quadrature 
to perform the calculation. Freericks et al noticed that 
averaging the results of two Gaussian quadratures with 
n and n + 1 points works better than choosing (2n + 1) 
points for the quadrature^ For the equilibrium case we 
adopt this trick. In Fig. [7] we plot the lesser Green 
function G<(T,t rel ) = J2k G< ( k \ T ^rei)/N calculated 
within the Kadanoff-Baym- Wagner NEDMFT in com- 
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FIG. 7: (Color online) The imaginary part of the equilibrium 
lesser Green function ImG < (T = 0,t re i) calculated within 
the Kadanoff-Baym- Wagner NEDMFT by using two Gaussian 
quadratures (GQ) of 50 and 51 points (black line) and by 
one Gaussian quadrature of 101 points (blue line) (At — 0.1, 
At — 0.1, to — —15) in comparison with the result calculated 
by solving the equilibrium DMFT equations in frequency (red 
line) . The inset focuses the imaginary part of the equilibrium 
lesser Green function G < (T = 0, t re i) obtained by performing 
the two Gaussian quadratures of 50 and 51 points in a small 
range of t re i for different At (At — 0.1, to = —15). The 
model parameters U = 0.5, j3 = 10. 
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parison with the one obtained by solving the equilibrium 
DMFT equations in frequency. Here we have converted 
the results from the time variables t and t' to Wigner 
average T = (t + t')/2 and relative t re \ = t — t' time vari- 
ables. The NEDMFT calculations are performed with 
two Gaussian quadratures with 50 and 51 points. We 
also plot the NEDMFT result which is obtained by per- 
forming only one Gaussian quadrature with 101 points. 
The results plotted in Fig. [7] confirm the notice of Fre- 
ericks et al. Indeed, for a range of small t re i the lesser 
Green function calculated within the NEDMFT fits per- 
fectly with the one obtained within the DMFT. However, 
for t re i nearby the time cutoffs, the NEDMFT results ex- 
hibit spurious features of a nodal form due to finite size 
effects of the numerical procedures. These spurious fea- 
tures are greatly reduced by employing the trick of two 
Gaussian quadratures. However, as we will see later, the 
spurious features do not appear in the nonequilibrium 
case where the electric field is finite. The lesser Green 
function obtained within the NEDMFT fulfils the sum 
rule very well, as shown in the inset of Fig [3 Indeed, 
ImG < (T, t re i = 0) w 0.5 in the weak interaction case. 
The sum rule of higher-order moments of the lesser Green 
function is also fulfilled because the lesser Green func- 
tion fits perfectly with the DMFT one nearby t re l = 0. 
However, if one numerically calculates the sum rule of 
higher-order moments, a deviation from the exact value 
may appear due to numerical derivative calculations from 
discretized points^ The possible deviation of the sum 
rule of higher order moment does not necessarily mean 
an inaccuracy of the Green function. It relates to the 
finite value of At which may be not small enough for 
performing numerical derivative calculations. 



In Fig. [8] we plot the imaginary part of the equilibrium 
lesser Green function calculated within the Kadanoff- 
Baym- Wagner NEDMFT and the equilibrium DMFT 
for U = 2. This value of U corresponds to the insu- 
lator phase. It shows that for strong interactions the 
NEDMFT results fit well with the ones of the equilib- 
rium DMFT for small t re i. The spectral sum rule of the 
lesser Green function is well fulfilled. Indeed, for U = 2 
ImG<(T = 0,t ret = 0) is equal to 0.5123 for At = 0.1, 
and is equal to 0.5021 for At — 0.05 in comparison with 
the exact value 0.5. Around the minima in the curve 
of the imaginary part of the lesser Green function small 
deviations appear. The deviations can be reduced by 
decreasing At. For strong interactions the numerical re- 
sults slightly deviate from the equilibrium values due to 
the finite discretization of the time variables. Neverthe- 
less, the NEDMFT calculations for the equilibrium case 
for both weak and strong interactions show the numerical 
techniques employed here are accurate and controllable. 

In the nonequilibrium case, when the electric field is 
finite, we notice that the use of two Gaussian quadra- 
tures for the integration with the double Gaussian den- 
sity of states gives almost the same result as the use of 
one Gaussian quadrature. In Fig. [5] we plot the imaginary 
part of the lesser Green function calculated by using two 
Gaussian quadratures of 50 and 51 points in comparison 
with the one calculated by using one Gaussian quadra- 
ture of 51 points for both weak and strong interactions. 
It shows that the results of both quadrature methods are 
almost identical. In contrast to the equilibrium case, in 
the nonequilibrium case there are not spurious features 
nearby the time cutoffs. We have also checked the re- 
sults with more points for the Gaussian quadrature (in 
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FIG. 8: (Color online) The imaginary part of the equilibrium 
lesser Green function ImG < (T = Q,t re {) calculated within the 
Kadanoff-Baym-Wagner NEDMFT for different At (At = 
0.1, to = —15). The two Gaussian quadratures with 50 and 
51 points are performed. The equilibrium DMFT calculation 
result is also presented (red line). The model parameters 
U = 2, f3 = 10. 



FIG. 9: (Color online) The imaginary part of the nonequi- 
librium lesser Green function ImG < (T = 0,t re i) calculated 
within the Kadanoff-Baym-Wagner NEDMFT by using two 
Gaussian quadratures of 50 and 51 points (solid lines), and 
by using one Gaussian quadrature of 51 points (symbols) for 
weak (U — 0.5) and strong (U — 2) interactions (At = 0.1, 
At = 0.1, t = -15, p=W,E= 1). 
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FIG. 10: (Color online) The imaginary part of the lesser 
Green function G < (T = 0,t re i) calculated within the 
Kadanoff-Baym- Wagner NEDMFT in the weak interaction 
case for various At. The inset plots the real part of the lesser 
Green function. The model parameters U = 0.5, E = 1, 
j3 — 10, to = — 15. n — 51 is used for the Gaussian quadra- 
ture and At = 0.1. 



particular, with n = 101), and with finer At (in partic- 
ular, with At = 0.05). It turns out that the numerical 
results are mostly sensitive to the real time discretiza- 
tion. In the following for numerical calculations we use 
the single Gaussian quadrature with 51 points and make 
the integrations over the imaginary time with At = 0.1 
for /3=10. 

In Fig. [TU] we present the lesser Green function in the 
weak interaction case for various real time discretizations 
At and a fixed to. It shows that the imaginary part of 
the lesser Green function quickly converges with At. It 
also indicates that the spectral sum rule of the lesser 
Green function is fulfilled well. Indeed, in Table I we 
list the value of the spectral sum rule of the lesser Green 
function for various values of At. In the Kadanoff-Baym- 
Wagner formalism the spectral sum rule which is numeri- 
cally calculated is fulfilled better than in the formalism of 
the contour-ordered Green function^ The spectral sum 



TABLE I: The spectral sum rule of the nonequilibrium lesser 
Green function at the average time T = calculated within 
the Kadanoff-Baym- Wagner NEDMFT for weak (U = 0.5) 
and strong (U = 2) interaction (E = 1). Various values of At 
are used (At = 0.1, to = -15, /? = 10). 
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At 
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exact 


sum rule 
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FIG. 11: (Color online) The imaginary part of the lesser 
Green function G K (T = 0, t re i) calculated within the 
Kadanoff-Baym- Wagner NEDMFT in the strong interaction 
case for various At. The inset plots the real part of the lesser 
Green function. The model parameters U = 2, E = 1, j3 = 10, 
to = —15. n = 51 for the Gaussian quadrature and At = 0.1. 



rule of the retarded and advanced Green functions is ful- 
filled well too. The real part of the lesser Green function 
converges with At is less quickly. However, it also con- 
verges well for small At. For weak interactions the nu- 
merical calculations solving the Kadanoff-Baym- Wagner 
NEDMFT equations work very well. For other physical 
quantities such as the electric current and the double oc- 
cupation their convergences with At are also good. In 
particular, we can obtain converged results at the limit 
At — > by using a Lagrange interpolation formula. 

In Fig. [11] we plot the lesser Green function in the 
strong interaction case for various real time discretiza- 
tions At and a fixed to . The imaginary part of the lesser 
Green function converges with At well. However, its 
spectral sum rule slightly deviates from the exact value, 
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FIG. 12: (Color online) The time dependence of the electric 
current I(t)/jo calculated within the Kadanoff-Baym- Wagner 
NEDMFT for various time cutoffs t m with fixed A = 0.1 and 
to = -5 (U = 2, E = 1 At = 0.1, = 10). 
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FIG. 14: (Color online) The time dependence of the con- 
tribution of the initial correlations to the electric current 
AP c (t)/jo for various E and U in the weak interaction case. 
The results are already scaled by a quadratic extrapolation 
with At = 0.1, 0.065, and 0.5 (to = -5, At = 0.1, = 10). 
The inset plots the scaled contribution of the initial correla- 
tions AP c (t)/jo for different initial times to (U = 0.5, E = 1, 
At = 0.1, = 10). 



FIG. 13: (Color online) The time dependence of the electric 
current I(t)/jo calculated within the Kadanoff-Baym- Wagner 
(KBW) NEDMFT (red line) in the weak interaction case for 
different electric fields. For comparison the result obtained 
from the Keldysh (K) NEDMFT (black line), and the con- 
tribution of the initial correlations to the current (green line) 
are also plotted. The data are already scaled with a quadratic 
extrapolation (At = 0.1, 0.065, and 0.05). The model param- 
eters U = 0.5, = 10, Ar = 0.1, and E = 1 (E = 2) for 
upper (lower) panel. 



as presented in Table I. The real part of the lesser Green 
function converges with At not so fast as in the weak 
interaction case. In general, for strong interactions the 
extrapolations of the numerical results of the Kadanoff- 
Baym- Wagner NEDMFT in the limit At — > require a 
more effort. Often in order to obtain reliable data we 
have to carry the numerical calculations with At smaller 
than the ones in the weak interaction case. 

One can notice that the numerical results of the 
Kadanoff-Baym- Wagner NEDMFT are independent on 
the maximal time t m when At and to are fixed. In- 
deed, we plot the electric current I(t) obtained from 
the Kadanoff-Baym- Wagner NEDMFT for various t m in 
Fig. 1121 For larger t m the time window is larger and we 
can observe the behaviors of the system at a more long 
time. But for larger t m the numerical calculations are 
also more expensive in time. We have to compromise 
the computation time and the need of the time window 
width. 

In Fig. [T3] we present the electric current calcu- 
lated within both the Kadanoff-Baym- Wagner and the 
Keldysh NEDMFT in the weak interaction case. In 
the Keldysh NEDMFT the initial correlations are ne- 
glected. The contribution of the initial correlations to 



the current is defined as the difference of the currents 
calculated within the Kadanoff-Baym- Wagner and the 
Keldysh NEDMFT. Certainly, the current has been pre- 
viously calculated within the NEDMFT of the contour- 
ordered Green function^ We find after extrapolating to 
At — > the results of the Kadanoff-Baym- Wagner for- 
malism agree well with the ones obtained within the 
contour-ordered Green function NEDMFT. This indi- 
cates the equivalence of the Kadanoff-Baym- Wagner and 
the contour-ordered Green function formalisms, as ex- 
pected. However, the Kadanoff-Baym- Wagner formalism 
represents the contour-ordered Green function in the ma- 
trix form, the elements of which are the physical Green 
functions. It is also similar to the Keldysh formalism. 
The spectral sum rule obtained within the Kadanoff- 
Baym- Wagner NEDMFT is fulfilled very well. The elec- 
tric current displays the Bloch oscillations, as noticed by 
Freericks et al£2- For small and large electric fields (for 
instance, E = 1) the current is monotonously damped to 
zero value. However, when the electric field increases fur- 
ther (for instance, E = 2) the current develops beats. As 
shown in Fig.[T31 in the weak interaction case the Keldysh 
and the Kadanoff-Baym- Wagner NEDMFT qualitatively 
give the same current. The contribution of the initial 
correlations to the current also oscillates with time in 
the same way as of the full current. When the current 
displays beats the initial correlation contribution displays 
beats too, as shown in Figs. Q1J] and [T3J In the weak in- 
teraction case the initial correlation contribution to the 
current is small in comparison with the full current. How- 
ever, it never vanishes except for t < when the electric 
field is absent and the current vanishes too. In the inset 
of Fig. [14] we plot the initial correlation contribution to 
the current for different initial times to- It shows that 
the results are independent on the initial time if it is 
far enough from t = 0. The initial correlation contribu- 
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FIG. 15: (Color online) The time dependence of the dou- 
ble occupation D(t) calculated within the Kadanoff-Baym- 
Wagner (KBW) and the Keldysh (K) NEDMFT in the weak 
interaction case for different initial times to- The results are 
already scaled by a quadratic extrapolation with At = 0.1, 
0.065, and 0.5. The dotted line is the double occupation in 
equilibrium (E = 0). The model parameters U = 0.5, E = 1, 
At = 0.1, 13 = 10. 

tion seems to be finite even when the initial time is in 
the remote past. However, the initial correlations do not 
qualitatively change the current properties in the weak 
interaction case. 

We also calculate the double occupation D{t) = 
(c^(t)c(t)f^(t)f(t)) which can be computed through the 
lesser Green function Qf(t,t) defined in Eq. (f6Tj) by 

D(t) = -in f Qf(t,t). (68) 

In Fig. [T5] we plot the time dependence of the double 
occupation D(t) calculated within the Kadanoff-Baym- 
Wagner and the Keldysh NEDMFT in the weak interac- 
tion case. Before the electric field is turned on (t < 0), 
the double occupation calculated within the Kadanoff- 
Baym- Wagner NEDMFT is constant in a good agree- 
ment with the equilibrium value obtained by solving the 
DMFT equations in frequency. However, the Keldysh 
NEDMFT results are quite different. Within the Keldysh 
NEDMFT, the double occupation starts from its non- 
interaction value at half filling (Dq = 0.25), and then 
decreases to a steady value. This steady value is little 
larger than the equilibrium value. At the initial time to 
the Keldysh formalism starts with noninteracting system 
and electron correlations are absent. The results show 
that before the electric field is turned on, the Keldysh 
formalism cannot restore full electron correlations of the 
system. The full electron correlations are essentially elab- 
orated from the initial correlations which come from the 
dynamics of the system in the imaginary time branch 
of the Kadanoff-Baym contour. The neglect of the ini- 
tial correlations also means the lack of the electron cor- 
relations even when the system is still in equilibrium. 
In Fig. [15] we also plot the double occupation for differ- 
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FIG. 16: (Color online) The time dependence of the dou- 
ble occupation D(t) calculated within the Kadanoff-Baym- 
Wagner (KBW) and the Keldysh (K) NEDMFT in the weak 
interaction case for different initial times to- The results are 
already scaled by a quadratic extrapolation with At = 0.1, 
0.065, and 0.5. The dotted line is the double occupation in 
equilibrium (E = 0). The model parameters U = 0.5, E — 2, 
At = 0.1, 0=10. 



ent initial times to. It shows that even when the initial 
time goes to the remote past, the lack of electron correla- 
tions still occurs in the Keldysh formalism. Only in the 
Kadanoff-Baym formalism when the initial correlations 
are taken into account, full electron correlations are ob- 
tained. After the electric field is turned on (t > 0) for 
weak and strong electric fields (for instance, E = 1) the 
double occupation first oscillates strongly, and then is 
damped into less regular oscillations. However, when the 
electric field increases further (for instance, E = 2) the 
double occupation regularly oscillates even at long time, 
as shown in Fig. [TBI This behavior is reminiscent to the 
beats of the electric current. In the weak interaction case 
the difference of the double occupations calculated within 
the Kadanoff-Baym- Wagner and the Keldysh formalisms 
is small in relation with their values. After the turning 
on of the electric field the difference becomes smaller. In 
the weak interaction case the Keldysh formalism qualita- 
tively describes the behavior of the double occupation. 

In equilibrium when the interaction U > \[2 the den- 
sity of states opens a gap at the Fermi energy, and the 
system is insulator. It distinguishes between the weak 
and strong interaction cases. In general, in the strong in- 
teraction case the numerical calculations slowly converge 
with At. Usually, we have to use more small values of 
At in order to get reliable results. In Fig. [17] we plot 
the electric current in the strong interaction case. In 
contrast to the weak interaction case, the current does 
not display the regular Bloch oscillations. The current 
oscillations are rather irregular and quenched. However, 
the current calculated within the Keldysh NEDMFT still 
exhibits the regular Bloch oscillations. This shows that 
the initial correlations are important in the strong in- 
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FIG. 17: (Color online) The time dependence of the electric 
current I(t)/jo calculated within the Kadanoff-Baym- Wagner 
(KBW) and the Keldysh (K) NEDMFT for strong interaction 
U = 2, E — 1. The contribution of the initial correlations is 
also plotted. The inset plots the initial correlation contri- 
bution for different initial times to- The results are already 
scaled by a cubic extrapolation with At = 0.05, 0.035, 0.025, 
and 0.02 ( At = 0.1, f3 = 10). 



FIG. 18: (Color online) The time dependence of the dou- 
ble occupation D(t) calculated within the Kadanoff-Baym- 
Wagner (KBW) and the Keldysh (K) NEDMFT in the strong 
interaction case for different initial times to- The results are 
already scaled by a cubic extrapolation with At = 0.05, 0.035, 
0.025, and 0.02. The dotted line is the double occupation in 
equilibrium (E = 0). The model parameters U = 2, E = 1, 
At = 0.1, 0=10. 



teraction case. They are a main factor for quenching 
the current oscillations. The contribution of the initial 
correlations to the current is not small as in the weak in- 
teraction case. It is of order of the current. In the inset 
of Fig.[T7]we also plot the initial correlation contribution 
to the current for different initial times to- It shows that 
the contribution remains the same as the initial time goes 
to the remote past. Thus, in the strong interaction case 
the initial correlations become significant, and dominate 
the overall properties of the current. The neglect of the 
initial correlations may cause artifacts in the nonequilib- 
rium properties of the current. 

In Fig.[TH]we plot the double occupation in the strong 
interaction case. Before the turning on of the elec- 
tric field, the time dependence of the double occupation 
within both the Kadanoff-Baym- Wagner and the Keldysh 
NEDMFT is similar to the weak interaction case. Within 
the Kadanoff-Baym- Wagner formalism the double occu- 
pation is constant for t < 0. The constant value is in 
a good agreement with the equilibrium value, although 
there are very little deviations due to the finite size effects 
in the numerical calculations. The double occupation 
calculated within the Keldysh formalism first starts with 
the noninteraction value Dq = 0.25 at the initial time to, 
then relaxes to a steady value. Like in the weak interac- 
tion case, the steady value is not the equilibrium value. 
It again indicates that the Keldysh formalism losses a 
some part of electron correlations. In the strong interac- 
tion case this lack of electron correlations becomes sig- 
nificant. As a consequence, after the turning on of the 
electric field the lack of electron correlations also remains 
significant. Due to the quenching of the Bloch oscilla- 
tions in the strong interaction case, the double occupa- 



tion reaches a steady value at a long time. The steady 
values obtained within the Kadanoff-Baym- Wagner and 
the Keldysh formalisms are quite different. They indicate 
the important contribution of the initial correlations. As 
shown in Fig. [18] the results do not change when the ini- 
tial time goes to the remote past. For strong interactions 
the Keldysh formalism losses a significant part of electron 
correlations both before and after the turning on of the 
electric field. It cannot correctly describe the nonequi- 
librium properties. 



VI. CONCLUSIONS 

In this paper we present the Kadanoff-Baym- Wagner 
formalism for nonequilibrium systems. The formalism 
is based on the Wagner representation of the contour- 
ordered Green function. Within the Kadanoff-Baym- 
Wagner formalism the Green functions satisfy the proper 
boundary conditions. The initial correlations essentially 
distinguish between the Kadanoff-Baym- Wagner and the 
Keldysh formalisms. We derive the system of equations 
for nonequilibrium Green functions, and solve it within 
the truncated and self-consistent perturbation theories as 
well as within the NEDMFT. As a benchmark we exam- 
ine the equilibrium FKM by the Kadanoff-Baym- Wagner 
NEDMFT. The results show a good agreement between 
the Kadanoff-Baym- Wagner NEDMFT in equilibrium 
and the equilibrium DMFT. In the nonequilibrium case 
the Green functions obtained within the Kadanoff-Baym- 
Wagner NEDMFT satisfy the spectral sum rule well. The 
derived Kadanoff-Baym- Wagner equations for nonequi- 
librium Green functions are an alternative useful method 
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for studying noncquilibrium systems. 

In this paper we also emphasize the initial correla- 
tions. Within the perturbation theory the initial cor- 
relations always finite even when the initial time goes 
to the remote past. The electric current calculated 
within the truncated perturbation theory shows that 
the Kadanoff-Baym- Wagner formalism overestimates the 
current, whereas the Keldysh formalism underestimates 
it. However, the Kadanoff-Baym- Wagner perturbation 
theory shows a better agreement with the exact solution. 
For a long time the truncated perturbation theory fails 
to describe the physical properties. The self-consistent 
perturbation theory gives better results than the trun- 
cated perturbation theory. The time domain in which the 
self-consistent perturbation theory gives reasonable re- 
sults is wider than the one in the truncated perturbation 
theory. However, the self-consistent perturbation theory 
cannot reproduce the beat behaviors of the current for 
very strong electric fields. Within the perturbation the- 
ories the initial correlations do not qualitatively change 
the perturbation results. Since both the Kadanoff-Baym- 
Wagner and the Keldysh perturbation theory results are 
close, use of the Keldysh approach is more convenient 
since its equations are simpler. In the infinite dimension 
limit the NEDMFT gives the exact solution. Examining 
the NEDMFT within both the Kadanoff-Baym- Wagner 
and the Keldysh formalisms one can figure out the role of 
the initial correlations. For weak interactions the initial 
correlations give only small contributions to the physical 
quantities such as the electric current or the double oc- 



cupation. However, they remain finite for the long time 
limit. The initial correlations are also important even 
before the electric field is turned on when the system is 
still in equilibrium. Without the initial correlations the 
system cannot restore the full electron correlations. For 
strong interactions the initial correlations become signif- 
icant, and dominate the physical properties. Without 
taking into account the initial correlations the Keldysh 
formalism can qualitatively describe the noncquilibrium 
properties of the system only for weak interactions. For 
strong interactions it fails to count full electron correla- 
tions. The neglect of the initial correlations may cause 
artifacts in the nonequilibrium properties of the system. 
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